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Expcrimcntal measurements of branching fractions 
and CP-violating charge asymmetries in rare B decays 
play an important role in testing the theoretical predic- 
tions of the standard model and its extensions. We re- 
port the results of branching fraction measurements for 
B° meson decays to r]K°, r]T], r]<p, -quo, rj'K , r/'r]', i]'(f>, 
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and rj'u final states and of branching fraction and charge 
asymmetry measurements for B + decays to r)ir + , r]K + , 
rj'n + , and rj'K + We search for charge asymmetry by 
measuring 



r- + r+ 



(i) 



where T ± = T(B ± — > f^) is the decay width for a 
given charged final state f^. These branching fraction 
and charge asymmetry measurements represent an im- 
provement over previous results published by BaBAR 
and Belle [3j. 

The branching fractions and charge asymmetries of 
the charmless hadronic B decays are predicted using ap- 
proaches based on QCD factorization [J, |a, |y, |7j and fla- 
vor SU(3) symmetry [1, H, U3]. These B decays proceed 
through loop (penguin) and suppressed tree diagram am- 
plitudes, as shown in Fig. [T] The branching fraction and 



B ' 




K ' 




(a) 



(b) 




u, c, t 




(c) 



(d) 



FIG. 1: Examples of Feynman diagrams involved in decays 
studied in this paper: (a,b) penguin diagrams, (c) Cabibbo- 
suppressed tree diagram, (d) gluonic penguin diagram. 



charge asymmetry measurements may provide sensitivity 
to the presence of heavy non-standard model particles 
in the loop diagrams 111] , The measured rf K branching 



fraction is found to be much larger than the r\K one 0,0]. 
Many suggestions have been proposed to explain such a 
difference, including flavor singlet enhancement [l2T |. in- 
trinsic charm [l3| , and constructively interfering internal 
penguin diagrams [3, EH- This last approach is sup- 
ported by next-to-leading order QCD factorization cal- 
culations [||. 

The CP- violating parameters S v >k and S^k, measured 
in the time-dependent analysis of rj'K and (j>K° de- 
cays [lH, are expected to equal S C c S ~ sin 2/3, where S C c S 
is measured in the Cabibbo-Kobayashi-Maskawa (CKM) 
favored b — > ccs decays, if penguin b — > s transitions 
are dominant. However, CKM-suppressed amplitudes 



and color-suppressed tree diagrams can introduce addi- 
tional weak phases whose contributions may not be neg- 
ligible @, m m. As 

a consequence, deviations from 
sin 2/3 may occur even within the standard model. Rates 
of the decay modes to 7777, r]4>, rfrf , and rf<j) are used in 
flavor SU(3)-based calculations of the \Scca — Sf \ (with 
/ = T]'K, 4>K) [l7| bound. This bound may be improved 
by more precise measurements of the branching fractions 
of these modes. 

The charge asymmetry is expected to be sizable in 
r]K + and suppressed in 77' K+ decays 0, ©, 03. 

How- 
ever, different approaches predict the two asymmetries 
to have the same @ or opposite [|[ signs; precise mea- 
surement of such asymmetries can discriminate between 
these models. Furthermore, the charge asymmetries in 
i]'tt + and r)ir + decays are expected to be sizable [y, 
with model-dependent predictions for their magnitudes. 

The results presented here are based on the full dataset 
collected with the BABAR detector [2(| at the PEP-II 
asymmetric-energy e + e~ collider located at the SLAC 
National Accelerator Laboratory. An integrated lumi- 
nosity of 426 fb _1 , corresponding to N B g = 467 x 10 6 
BB pairs, was recorded at the T(AS) resonance (center- 
of-mass energy ^fs = 10.58 GeV). A further 44 AT 1 
was collected approximately 40 MeV below the resonance 
(off-peak) for the study of the e + e~ — > qq background, 
where q is a u, d, s, or c quark. 

Charged particles are detected, and their momenta 
measured, by a combination of a vertex tracker, consist- 
ing of five layers of double-sided silicon microstrip detec- 
tors, and a 40-layer drift chamber, both operating in the 
1.5 T magnetic field of a superconducting solenoid. We 
identify photons and electrons using a CsI(Tl) electro- 
magnetic calorimeter (EMC). Further charged-particle 
identification (PID) is provided by the average energy 
loss (dE/dx) measurements in the tracking devices and 
by the information provided by an internally reflecting 
ring-imaging Cherenkov detector (DIRC) covering the 
central region. 

K®, w, and 



P , "SI 

1 77 {Vm 



7T 
7T+7T 



candidates 



7T 



A 



We select 77, 77', 4 
through the decays 77 
77' -> tjtt+tt- with rj -> 77 (77^), V 

4> -> K+K~, p° -> TT+TT", K° -> 7T+7T-, W -> 7T+ IT , 

and 7T° — > 77. We do not study the decay B° — > 7/77' 
with both 77' mesons decaying to fry, because it suffers 
large backgrounds. Requirements applied to the photon 
energy E-y and to the invariant mass of the B daughters 
are listed in Table [H The requirements on the 77 and 77' 
invariant masses depend on the decay mode. Branching 
fractions of charged B decays with 77' in the final state 
and of B° — > r/'K® are higher than those of the other neu- 
tral B modes. In neutral decay modes we apply a tighter 
requirement on 773^- invariant mass in order to prevent 
possible contamination from BB background. The dif- 
ferent requirements on the 77' mass increase the purity of 
the charged B and 77' K® modes and enhance the selec- 
tion efficiency for the other neutral B decay modes. The 
energy (momentum) of the 7T° (77) candidates is required 
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TABLE I: Selection requirements on the invariant masses of 
the signal resonances and on the laboratory energies of the 
photons coming from their decays. 



State 


Invariant mass 


(MeV/C 3 ) E 


y (MeV) 


n 

7T 


120 < 




< 150 


> 30 


Prompt 77 77 


505 < 




< 585 


> 100 


Secondary r? 77 


490 < 


TH-y-y 


< 600 


> 50 f 


773tt in B + decays 


534 < 




< 561 


— 


rjz-n in B° decays 


535 < 




< 555 




rj^^ in B + and 


945 < 




< 970 




B° ri'K° s decays 










ri' vlT7r in other B° decays 


930 < 




< 990 




rj' pl in B + decays 


930 < 




< 980 


> 200 


rf„ in B° - rf„K% 


930 < 


TTip-y 


< 980 


> 100 


r)' pl in other B° decays 


910 < 


TfXp^ 


< 990 


> 200 




470 < 




< 990 




U) 


735 < 




< 825 






1012 < 


m K+K- 


- < 1026 




K° s 


486 < 


TTIti-ti- 


< 510 





f £ 7 > 100 MeV in the B + -> T]' vnv K + and B + -> t/^tt 
decay modes. 



to exceed 200 MeV (200 MeV/c) in the laboratory frame. 
The prompt charged tracks in B + — > ?/7r + and secondary 
charged tracks in 77, r/' , and u> candidates are required to 
have DIRC, dE/dx, and EMC signatures consistent with 
the pion hypothesis. After selection, we constrain the 77, 
77', and 7T° masses to their world average values [2l|. The 
prompt charged track in B + — > rj K + is required to be 
consistent with the kaon hypothesis. The signatures for 
the charged kaons from <fr decays are required to be incon- 
sistent with hypotheses for electrons, pions and protons. 
For the prompt charged track in B + decays to r)K + and 
7/7r + , we define the variables Ck and C n as 

ameas _ a ex P 
~ „meas ' \ Z ) 

where 9%™ s (6%*) is the measured (expected) DIRC 
Cherenkov angle and cr^ e ^ s is its uncertainty, for the 
kaon and pion hypothesis, respectively. We require 
-3 < C K < 13 and -13 < C w < 3. For K° s candidates 
we require a vertex \ 2 probability larger than 0.001 and 
a reconstructed decay length greater than three times its 
uncertainty. 

We reconstruct the B meson candidate by combining 
the four-momenta of the final state particles and impos- 
ing a vertex constraint. A B meson candidate is kine- 
matically characterized by the energy-substituted mass 
rriES = \J s/4 — p% and energy difference AE = Eb — 
\\fs, where (Eb,Pb) is the S-meson four-momentum 
vector expressed in the T(45) rest frame. For sig- 
nal events, the Toes an d AE distributions peak around 
5.28 GeV/c 2 and zero, respectively. We require 5.25 < 
toes < 5.29GeV/c 2 and \AE\ < 0.2 GeV for all decay 



modes except B° — > r/K®, where we require —0.15 < 
AE < 0.2 GeV in order to suppress most of the back- 
ground from radiative B decays. 

Backgrounds arise primarily from random combina- 
tions of tracks and neutral clusters in e + e _ — > qq contin- 
uum events. We use large samples of Monte Carlo (MC) 
simulated [22| events and control samples to optimize cri- 
teria to suppress the background. We reject continuum 
events by using the angle #t between the thrust axis of 
the B candidate in the T(AS) frame and that of the rest of 
the event. The thrust axis of the B candidate is given by 
the thrust axis of the B decay products. The distribution 
of I cos6*t| is sharply peaked near 1.0 for jet-like qq pair 
events and is nearly uniform for T(AS) — > BB events. We 
require | cos£?t| < 0.9 (0.85 for f]' pl ^ + ', 0.8 for 7y 77 w and 
r]'u>). To discriminate against r-pair and two-photon 
backgrounds, and to better describe the event shape, we 
require the event to contain at least three charged tracks, 
or one track beyond the minimum required for the signal 
decay topology, whichever is larger. 

In r] — > 77 (0) decays, we define Tt v (H<f>) as the cosine 
of the angle between the direction of a daughter 7 (K) 
and the flight direction of the parent of 77 ((j)) in the 77 (<p) 
rest frame; for 77' , TL P is the cosine of the angle between 
the direction of a daughter pion and the flight direction 
of the 77' in the p rest frame. For B decays containing an 
uj meson in the final state we define H u as the cosine of 
the angle between the B recoil direction and the normal 
to the plane defined by the uj daughters in the uj rest 
frame. We require \H V \ < 0.95 in B° — > 7777 decay modes. 
We reject candidate events if \H P \ > 0.9 (> 0.75 in the 
B + — > r] / pi n + decay mode). 

For the B° -> t] 17 K q s (B+ -> ?7 7 7 ^+, h+ = K + ,n+) 
decay, the main source of BB background is the B° — » 
tt°K® (B + — * ir°h + ) decay. To suppress this back- 
ground, we search for ir° candidates with a photon in 
common (overlapping) with the 77 candidate from the re- 
constructed signal B candidate. We require the ir° mass 
not to be in the range (0.117, 0.152) GeV/c 2 for the 
B° -> ri 17 K° decay mode, and (0.118, 0.150) GeV/c 2 
for the B + — > 77 77 /i + decay modes. Further suppression 
of this background is obtained with suitable requirements 
on \H n \ and on the energy of the second (non-overlapping 
with 77) 7T° photon (E^ nd ). We optimize these require- 
ments by maximizing S/VS + B, where S (B) is the 
number of signal (background) events surviving the selec- 
tion. We find the optimal criteria to be \TC V \ < 0.966 and 
E^ nd < 0.207 GeV for the B° -> 77^° decay mode, and 
\H V \ < 0.977 and E^ nd < 0.143 GeV for the B+ ~* r,^h+ 
decay modes. 

We find a mean number of B candidates per event in 
the range 1.0-1.4, depending on the final state. Signal 
events are divided into two categories: correctly recon- 
structed (CR) signal where all candidate particles come 
from the correct signal B, and self cross-feed (SCF) signal 
where at least one candidate particle is exchanged with 
a particle coming from the rest of the event. Simula- 
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tions show that the fraction of SCF candidates is in the 
range (3-7)% in charged B decay modes and (2-20)% 
in neutral B decay modes. If an event has multiple B 
candidates, we select the candidate with the highest B 
vertex \ 2 probability, determined from a vertex fit that 
includes both charged and neutral particles 23]. This 
algorithm selects the correct candidate, if present, with 
an efficiency of (91-99)% and introduces negligible bias. 

We obtain yields from unbinned extended maximum- 
likelihood (ML) fits. The main input observables are AE, 
i«es, and a Fisher discriminant T (24j. Where relevant, 
the invariant masses m res of the intermediate resonances 
and angular variables Ti. are used. The Fisher discrimi- 
nant T combines five variables: the angles with respect to 
the beam axis of the B momentum and B thrust axis, the 
zeroth and second angular moments Lo,2 of the energy 
flow about the B thrust axis, and the absolute value of 
the continuous output of a flavor-tagging algorithm. The 
first four variables are evaluated in the T(4S) rest frame. 
The moments are defined by L r = J2 S Ps x l cos @s | r , where 
6 S is the angle with respect to the B thrust axis of track 
or neutral cluster s with momentum p s , and the sum 
excludes the B candidate. Flavor tagging information is 
derived from an analysis of the decay products of the non- 
signal candidate B meson (Stag), using a neural network 
based technique [25[. The output value of the tagging 
algorithm reflects the different final states identified in 
Stag decay. In particular, the presence of a lepton in the 
final state usually results in a large tagging output value, 
for both B°B° and B + B~ events. Since leptons are not 
generally present in continuum background events, the 
inclusion of the tagging algorithm output in T improves 
its discriminating power between continuum background 
and BB events. The coefficients of T are chosen to max- 
imize the separation between the signal and the contin- 
uum background. They are determined from studies of 
MC signal events and off-peak data. 

The set of probability density functions (PDF) used in 
the ML fits, specific to each decay mode, is determined 
on the basis of studies with MC samples. We estimate 
BB backgrounds using MC samples of B decays. Where 
needed, we add components to account for BB back- 
ground events with a Toes or AE distribution that peaks 
in the signal region and for background from B meson 
decays with charmed particles in the final state. 

The extended likelihood function is 



N 



''{Pj ( X «) 



(3) 



where N is the number of input events, rij is the number 
of events for hypothesis j (j = 1 for signal, j = 2 for 
continuum background, and j ' = 3 for BB background), 
and Vj(xi) is the corresponding PDF evaluated with the 
observables Xj of the i th event. In the B° — ► rj'ui, rfcf), and 
il' pj uj decay modes the signal includes both the CR and 
SCF signal components with the SCF fraction fixed to 
the value estimated from simulation. Due to the similar 



kinematics and branching fractions of the r/K + and tot + 
decay modes, we perform a combined fit to extract the 
two signal yields and charge asymmetries. In this fit we 
use the Ck and variables to discriminate the mass 
hypothesis of the prompt track. Since the correlations 
among the observables in the data are small, we assume 
each Vj to be the product of the PDFs for the separate 
variables. Correlations between the rjK + and tot + signal 
yields (charge asymmetries) are below 5% (7%). 

We determine the PDF functional form and parame- 
ters from MC simulation for the signal and BB back- 
grounds, and from sideband data (5.25 < toes < 5.27 
GeV/c 2 ; 0.1 < |AS| < 0.2 GeV) for the continuum back- 
ground. For B + — > rjh + decay modes, PDF functional 
form and parameters for the continuum background are 
determined using off-peak data. We parameterize each of 
the functions Pi (toes), 'Pi(AE), Vj{!F), and the peak- 
ing components of Vj(m res ) with either a symmetric or 
a bifurcated Gaussian, the sum of two symmetric or bi- 
furcated Gaussian shapes, a bifurcated Gaussian distri- 
bution with exponential tails [26| or a Crystal Ball func- 
tion [27]], as required to describe the distribution. Slowly 
varying distributions (m rcs and AE for the continuum 
background, and angular variables) are represented by 
linear or quadratic functions. For the continuum back- 
ground, the Toes distribution is described by the ARGUS 
function [28[ . Large data control samples of B decays to 
charmed final states of similar topologies are used to ver- 
ify the simulated resolutions in toes and AE. Where 
the control samples reveal differences between data and 
MC samples in mass (energy) resolution, we correct the 
mean and scale the width of the mass (energy) distribu- 
tion used in the likelihood fits. 

The validity of the fit procedure and PDF parameter- 
ization, including the effects of unmodeled correlations 
among observables, is checked with simulated experi- 
ments. This is done by embedding a number of signal 
and peaking BB background events from fully simulated 
MC samples and by drawing a number of qq and charm 
BB events from PDFs, according to the values found in 
the data. In each fit the free parameters are: the yields, 
the charge asymmetry for the signal and continuum back- 
ground, and several parameters describing the AE, TOes> 
and T distributions of the continuum background. A sys- 
tematic uncertainty due to fixing signal and background 
parameters in the fit is accounted. The charge asym- 
metry for BB background is fixed to zero in the fit. A 
systematic is evaluated to account for this restriction. 

Table ED and Table DID show, for B° and B+ decays, 
respectively, the measured yields, fit biases, efficiencies, 
and products of daughter branching fractions for each 
decay mode. The efficiency is calculated as the ratio of 
the number of signal MC events after the event selection 
to the total generated, and is corrected for known dif- 
ferences between simulations and data. We compute the 
branching fractions from the fitted signal event yields, 
reconstruction efficiencies, daughter branching fractions, 
and the number of produced B mesons N B g, assum- 
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TABLE II: Fitted signal event yield and fit bias in events (ev), detection efficiency e, daughter branching fraction product 
Y[Bi, significance <S, and measured branching fraction B with statistical error for each B° decay mode. For the combined 
measurements we give the significance (with systematic uncertainties included) and the branching fraction with the statistical 
and systematic uncertainties (in parentheses the 90% CL upper limit). Significances greater than 7 standard deviations (a) are 
omitted. 
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Yield (ev) 


Fit bias (ev) 


€ (%) 


n ft (%) 
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1.4 
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1.8 


0.59+SiI 

0.40 
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490^24 


-2 
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6.1 




64.9+3^ 
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28.3 
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72.4 ± 3.0 














68.5 ± 2.2 ± 3.1 
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2.2 
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VrjiYir np~y 
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0.6t° : ? 
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r\ 'ui 
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+1 
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3.4 
1.2 
3.6 


1 no+u.04 
l.UO_ 46 

0.94^-j 1 
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TABLE III: Fitted signal event yield and fit bias, detection efficiency e, daughter branching fraction product nft> measured 
branching fraction B, charge asymmetry Ach with statistical error, and significance Sa of the charge asymmetry for each charged 
decay mode. For the combined measurements we give the branching fraction, the charge asymmetry and the significance of the 
charge asymmetry with the statistical and systematic uncertainties. 



Mode 


Yield (ev) 


Fit bias (ev) 


e(%) 


rift (%) 


6(10~ 6 ) 


Ach 


Sa(ct) 


?7 77 7T + 


286t 3 3 i 
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0.3 
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3.3 






+1 
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-0.25 + 0.19 


1.3 




llltg 


+7 
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z - a -0.8 
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2.1 












3.5 ± 0.6 ± 0.2 


+0.03 ± 0.17 ± 0.02 


0.2 




1601+11 


-5 


28.7 


17.5 


68.5li;« 


-0.004 + 0.027 


0.2 


29911^ 


-10 


29.3 


29.4 


74.6 + 1.8 


+0.016 + 0.023 


0.7 












71.5 ± 1.3 ± 3.2 


+0.008lo;"g ± 0.009 


0.4 



ing equal production rates of charged and neutral B 
pairs from T(4S) decays. We correct the yields for any 
bias measured with the simulations. We combine re- 
sults from different sub-decays by adding the values of 
-2 In (£/£ 

max) (parameterized in terms of the branch- 
ing fraction or charge asymmetry), where C max is the 



value of £ at its maximum, taking into account the cor- 
related and uncorrelated systematic errors. We report 
the branching fractions for the individual decay chan- 
nels and their significances S in units of standard devia- 
tions (a). For B° — > r]'K® and all charged decay modes, 
where the significance of the branching fraction is always 
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greater than 7c, the value of S is omitted. For the com- 
bined measurements we also report the 90% confidence 
level (CL) upper limits of the branching fraction for the 
B° modes where the significance is less than 5cr. For 
charged B decays we give the combined result for the 
charge asymmetry A c h and its significance 5a in units of 
a. 

The statistical uncertainty on the signal yield and 
charge asymmetry is calculated as the change in the cen- 
tral value when the quantity — 21n£ increases by one 
from its minimum. The significance is calculated as the 
square root of — 2 In (Co/Cmax), with systematic uncer- 
tainties included, where Cq is the value of L for zero signal 
events or zero value for the charge asymmetry. We deter- 
mine a Bayesian 90% CL upper limit on the branching 
fraction, assuming a uniform prior probability distribu- 
tion, by finding the branching fraction below which lies 
90% of the total of the likelihood integral in the positive 
branching fraction region. 

Figures [2] and [3] show the projections onto the tties and 
AE variables for the four neutral decay modes that have 
a branching fraction significance greater than 3d, and for 
the four charged decay modes, respectively. For each de- 
cay mode we optimize a requirement on the probability 
ratio Vi/(Vi+V2+'P3) in order to enhance the visibility 
of the signal. The probabilities Vj are evaluated with- 
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FIG. 2: The B° candidate ities and AE projections for r/Ks 
(a, b), rjuj (c, d), rj'Kg (e, f), and rfu) (g, h) decays, with sub- 
decays combined. Points with errors represent the data, solid 
curves the full fit functions, and dashed curves the background 
functions. 
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FIG. 3: The B + candidate tties and AE projections for 
?77r + (a, b), r]K + (c, d), r]'-K + (e, f), and rj' K + (g, h) decays, 
with sub-decays combined. Points with errors represent the 
data, solid curves the full fit functions, and dashed curves the 
background functions. 



out using the variable shown. The points show the data 
that satisfy such a requirement, while the solid curves 
show the total rescalcd fit functions. In ifuj decays, a 
fit performed on uj mass sidebands TO W7r7r < 735 MeV/c 2 
or m TWlr > 825MeV/c 2 shows that contamination from 
possible B° — > 77'7r + 7r~7r° background is negligible. 

The main sources of systematic error include ML fit 
bias (0-14 events) and uncertainties in the PDF param- 
eterization (0-12 events). The ML fit bias systematic 
error is taken to be half of the bias, summed in quadra- 
ture with its statistical uncertainty. The uncertainties re- 
lated to the PDF parameterization are obtained by vary- 
ing the PDF parameters within their errors. Published 
world averages [2l| provide the uncertainties of the B- 
daughter branching fractions (0-4)%. These uncertain- 
ties are the main contribution to the systematic errors of 
the B — > rf K decay modes. The uncertainty on N B g is 
1.1%. Other sources of systematic uncertainty are track 
(1%) and neutral particle (3-6%) reconstruction efficien- 
cies; selection efficiency uncertainties are 1% each for the 
cos#t and PID requirements. Using large inclusive kaon 
and B decay samples we estimate a systematic uncer- 
tainty for A c h of 0.005 due to the dependence of the 
reconstruction efficiency on the charge of the high mo- 
mentum K . Other sources of systematic uncertainties 
for Ach are the fit bias (0-0.02) and the presence of a fit 
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bias in the signal yield (0.02-0.03). The systematic un- 
certainty due to fixing the value of the charge asymmetry 
in BB background components is taken to be the largest 
deviation observed when varying this value of ±10%, and 
is in range (0-0.02). 

In summary we present updated measurements of 
branching fractions for eight B° and four B + decays to 
charmless meson pairs. The results shown in Table [TT1 and 
Table HIT] are consistent with, but generally more precise 
than, previous measurements @, y[ and supersede our 
previous ones Q. The branching fraction results are in 
agreement with predictions within the theoretical uncer- 
tainties that limit discrimination between different mod- 
els 0, i, i, 0, i, ©, 03. We find evidence for three B° 
decay modes: rjK° (3.5c), tjuj (3.7a) and r)'u> (3.6a). In 
the decay mode B + — > r]K + we find evidence at 3.3<r for 
non-zero charge asymmetry, in agreement with theoret- 
ical predictions E^l- Discrimination between QCD 
factorization [|| and flavor SU(3) [9] symmetry models, 
based on the relative sign of the charge asymmetry in 
B + — > r)K + and B + — ► r)'K + decays, is limited by the 
accuracy of the latter measurement. The measurement 
of A c h for rj'ir + shows a slightly better agreement with 
the QCD factorization prediction @ than with the flavor 
SU(3) symmetry based model 0], within large theoretical 



and experimental uncertainties. 
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